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Abstract 

In the present paper, we have briefly reviewed Sznajd's sociophysics model and its 
variants, and also we have proposed a simple Sznajd like sociophysics model based on 
Ising spin system in order to explain the time evaluation of resistance probability of a 
closed community against occupation. Using a numerical method, we have shown that 
time evaluation of resistance probability of community has a non-exponential character 
which decays as stretched exponential independent of the number of soldiers in a one- 
dimensional model. Furthermore, it has been astonishingly found that our simple 
sociophysics model belongs to the same universality class with random walk process on 
the trapping space. 
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1 Introduction 



Modelling of the some social phenomenon using the Ising spin system has been of 
considerable interest in the statistical physics. One of the important study in this 
area is undoubtedly the Sznajd model pQ which deals with opinion evolution in closed 
community. In the Sznajd model, each site of a one- or two-dimensional lattice carries 
a spin that can be either up (Republican) or down (Democrat) and represent one of 
two possibilities on any questions. Two neighboring parallel spins i.e., two neighboring 
people sharing the same opinion, convince their neighbors of this opinion. If they do not 
have the same opinion, then either they do not influence their neighbors or they bring 
their neighbors to the opposite position. In this model, the system evolves from one 
time step to another through a random sequential updating mechanism and it always 
reaches an overall consensus, i.e., the system ends up in a fixed point: either all spins 
point up or all point down. 

Many variants of this Sznajd model are studied. Various rules are defined and 
discussed in [2j. In all of them, spins pairs are selected randomly for trying to convince 
their neighbors. Stauffer stated that most variants of the systems always end up in 
a fixed point, where all spin are the same, or where the opinions are ordered anti- 
ferromagnetically in one or two dimensional for these rules (Sj. 

Memory of the system had been discussed based on the Sznajd model with syn- 
chronous updating in Ref. 4 . In the case of a system corresponding to a community 
always going to an overall consensus, this situation causes a dictatorship (all spins par- 
allel) in the Sznajd model. Therefore, to avoid dictatorship one can introduce a small 
number of non-conformist which are not convinced by the Sznajd rule (Schneider |5|), 
or the Sznajd model can be studied in disordered lattices so that this is more realistic 
since human society does not follow a square lattice with every person having exactly 
four neighbors. The simplest case is a diluted square lattice, where every site is either 
empty or having a single spin To simulate political election results, the Sznajd 
model was generalized from two different opinion to number q of opinions 6 and also 
Stauffer discussed Sznajd model with number q of opinions. He stated that there is 
reachable consensus for q < 3 but not for q > 4 7 . Furthermore, a different persua- 
sion process with a continuous spectrum of opinions between and 1 were studied by 
Deffuant et al fSj. 

Another generalization of Sznajd model to a triangular lattice with spreading of 
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mixed opinion and with the pure antiferromagnetic opinion was studied by Chang 
who found the fixed points for all values of initial concentration of down spins 
for the mixed case of ferromagnetic and antiferromagnetic opinions. Furthermore, a d- 
dimensional Sznajd model of consensus finding process was simulated with 1 < d < 4 for 
the two possible opinions, and the density of never changed opinion (which correspond 
to persistent spins in Ising model |ll)j ) during the Sznajd consensus-finding process 
decays with t as t 1 / 9 were discussed. It was found that exponent 9 is compatible with 
the Ising value for one dimension although it is not compatible for higher dimensions 

en 

Although, political cases that have two state opinion are quite extreme example, 
such phenomena appears when a referendum has been made in a community or in a 
country. People vote either "Yes" or "No" in a referendum such as election either 
Democrats or Republicans in US political conditions. 

Consequently, it seems that the Sznajd model T and its variants [211211110101111121 
ITT71 ITT] suggests to be able to make the physical model of some sociological or political 
events like Galam's sociophysics models [12| 113] , although there were resistance against 
modelling of the some sociophysics events from physicists [13] . 

The intent of this study motivated from Sznajd's 1 , Stauffer's jjllj and Galam's 
d8j papers is to introduce a simple sociophysics model to explain the time evaluation 
of resistance probability of a closed community against occupation. 



2 The Model 

Everybody knows that human history has been filled with conflicts, invasion or wars 
among countries or community. In historical age, either some countries (or community) 
have taken over other countries or some countries (or community) have been occupied 
by other countries in a military or different way. It is difficult to show that a country 
or a society had been able to stay out of such a process. 

Of course, conflicts, invasions, occupations and wars are quite complicate phenomen 
which generally occur on the economical, military, sociological, political, humanity 
levels etc. Therefore, it is impossible to suppose a sociophysics model comprised of all 
events in a conflicts, an invasion or a war. 

Nevertheless, in order to explain the dynamics of some socio-political events, a 
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sociophysical model can be set based on Ising spin systems just like Sznajd model even 
though it has to contain some definite limitations. In this study is considered a process 
in which a community (or a country) is under military occupation. To modelling such 
an occupation processes, some assumptions, of which details are given as follows, are 
considered. 

Initially, let us consider one-dimensional discrete lattice. Each site carriers a spin, 
S, which is either up or down as shown in Fig^ Host people of community and soldiers 
are represented with spin up (+1) and spin down (-1) respectively. 

+ 1 -1 +1 -1 +1 +1 +1 +1 -1 +1 +1 +1 



Fig. 1: Schematic representation of occupation process. Host people and occupying 
soldiers are represented with spin up (+1) and spin down (-1) respectively. Soldiers 
are distributed randomly among host people in the one-dimensional discrete space. 

Down spins are distributed randomly on the lattice governed by volume exclusion 
principle as seen in Fig^ This means that only one person (host or soldier) occupies 
a site. The number of soldiers is only a small fraction of the total population and then 
the density of soldiers, p (for one-dimensional model), is defined as number of soldiers 
over total population. 

We assumed that there is a over all consensus among host people in the community 
against occupation even if some exceptions exist. One expects that host people obey 
to this consensus at least initially. In this sense, community behaves as polarized 
at zero social temperature |15| against occupation just like Ising ferromagnet at zero 
temperature. 

As known, every person may be influenced by others, since they interact with each 
other in the community. Change of opinions of people can be simply explained by social 
impact [TH] . 

But in our model, soldiers have not been influenced by the host people. Their 
opinion about justifying the occupation does not change during the occupation process, 
since they may be stubborn, stable or professional etc. In this sense, it may be assumed 
that soldiers behaves like persistent spins in Ising spin system ^U], which are change 
their orientations corresponding to their opinions. In addition, mission of the soldiers 
is not only occupation, at the same time, they also want to convince the host people 
about justification of occupation. 
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Unlike soldiers, it can be conjecture that host people are influenced by soldiers even 
though they against occupation. They are exposed to intensive biased information or 
propagation, which play a role as persuasion mechanism. 

To explain the time evolution of one-dimensional sociophysics models, let us con- 
centrate on how our mechanism progresses. At this point, we shall define a resistance 
probability to occupation for all people (host or soldier) at any site i. For any soldier, 
the resistance probability, W%, could be zero during an occupation process since they 
believe in the justification of occupation. On the other hand, for host people at any site 
i, the resistance probability Wi is equal to one initially since they obey the community 
consensus. 

In this manner, soldiers affect neighbors who are host people, in other words, a 
host at any site i is influenced by the neighbors. Effected people may change their own 
opinions depending on resistance probability of the nearest neighbors about occupation. 
Namely, the phrase of 'winner takes all' is invalid herein, they might take only a 
probabilistic manner. Such a mechanism will depolarize the polarization (resistance 
probability) of all host people. 

For the any time step, resistance probability W% of any host people at site i is 
determined by probability of the nearest neighbors at the previous time as 



Eq.(^Q) clearly indicate that past plays an important role in this model. Also, the system 
evolves from one time step to another through synchronized updating mechanism. 

Total resistance probability of community will be depolarized by soldiers even if 
we assumed that community has a consensus about resistance to occupation initially. 
Total polarization for such a system could be calculated over every host people for any 
time step N as 



where mo is the initial number of host people, r represents any configuration and the 
sum is over all sites. If Eq.(j2J is averaged over different soldier configurations, then it 
can be rewritten as 




(1) 




(2) 




(3) 



r=l 



where R is the number of different configurations. 
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3 Numerical Results 



As an example, the schematic representation of the time evolution of resistance prob- 
ability due to Eq.Q is written using periodic condition for every spin (host people 
and soldiers) is shown in FigEJ In this figure, the first line represents the occupied 
community, and second line corresponds to its resistance probability at time N = 0. 
While probability values for soldiers remain unchanged, it decreases for host people as 
seen in Fig[2J 
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Fig. 2: Schematic representation of the time evolution of the resistance probability for 
every host people in the one dimensional occupied community model. 

Evolution of resistance probabilities for each host in time N is partly illustrated 
in FigEl Fortunately, Eq.® allow us to work with large lattice and over many differ- 
ent configurations. For this reason, Eq.Q has been solved numerically for lattice size 
L = 10.000, and density of the soldier p = 0.1, number of the time step iV = 5000 and 
number of the independent configuration R = 1000 with periodic boundary configura- 
tion by using a simple Fortran program. 

In order to understand clearly the time evolution of one-dimensional sociophysics 
model under conditions mentioned above, the numerical data were plotted as In < 
P(N) > versus N in FigEk, as log < P(N) > versus logiV in FigEb, as < P(N) > 
versus In N in FigOfc, and also as log(— In < P(N >) versus log N in FigOJl. 

It is explicitly seen from FigEK, FigOh) an d FigEt that there are no exponential, 
power, logarithmic laws dependence in the numerical data, respectively. But, as seen 
in Fig|21i, numerical data well fit to stretched exponential function as 

< p (N) >= Ae~ aNl3 (4) 
where A normalization factor, a relaxation constant, and (3 is the decay exponent of 
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Fig. 3: Numerical data for p = 0.1 were plotted In < P(N) > versus N in Fig.3a, 
log < P(N) > versus log TV in Fig.3b, < P(N) > versus IniV in Fig.3c, and log(— In < 
P(N) >) versus logA^ in Fig.3d. Solid-dot lines indicate numerical data and solid lines 
represent fitting curves in all figures. 
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the averaged resistance probability. 

In addition, it should be tested whether FigEJl satisfies to stretched exponential 
or not. If prefactor A is equal to 1 in Eq.Q the Fig02i would work as stretched 
exponential. However, if A is less than 1, small prefactor may give the impression of 
stretched exponential form, even for (3=1. Therefore, it can be plotted — In < P(N) > 
versus suitable powers of N, like iV 1 / 2 , iV 1 / 3 , etc., and find out the best straight line 
among the powers of N for long times (email from D. Stauffer). Hence, — In < P(N) > 
was plotted versus powers of N for p = 0.1 then the best straight fitting line for long 
times was obtained for (3 = 0.4 and (3 = 0.5 as seen in Fig0^ and FigHJ), respectively. 
These results confirm to this method used to find out stretched exponential exponents 
in FigEJi. Also, this test indicates that prefactor in Eq.Q does not effect results 
presented in this paper. 




Fig. 4: (a) — In < P(N) > versus iV ' 4 for p = 0.1, and (b) — In < P(N) > versus 

AT0.5 

for p = 0.1. Solid-dot lines indicate numerical data and solid lines represent fitting 
curves in all figures. 

To see whether the exponent /3 depends on p or not, we plotted (3 versus p in Fig03 
This figure shows clearly that exponent (3 is a function of density of soldiers. As seen 
from Fig|SJ (3 almost linearly decreases. It is inferred from Fig|3 that the averaged 
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Fig. 5: Change of the exponent (3 depending on density of soldiers p for one-dimensional 
model. Solid-dot line and dot line represent numerical data and fitting curve respec- 
tively. 

resistance probability < P{N) > has stretched exponential character independent of 
p. But, one cannot guaranty that it is true for all p values. As seen in Fig|Sl while j3 
monotonically decreases for small p whereas oscillating occurs in (3 for large p values, 
which indicates that fatal statistical errors may occur for large p values. 

4 Conclusions 

In all reasoning, it is inferred that numerical results obtained from simple sociophysics 
model studied herein obey a stretched exponential law independent of the number of 
soldiers. Also, we find out exponent depend on number of soldiers. These are quite 
remarkable aspects. 

Stretched exponential behavior indicates mathematically that decay for the rela- 
tively short times is fast, but for relatively long times it is slower. One can observe that 
this mathematical behavior corresponds to occupation processes in the real world. In 
generally, a military occupation is realized after a hot war. The community does not 
react to occupation since it occurs as a result of defeat. People are affected easily by 
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propaganda or other similar ways. Therefore, it is no surprise that resistance probabil- 
ity decrease rapidly at relatively short times. On the other hand, spontaneous reaction 
may begin against occupation in the community after the shock. Hence, community 
begins by regaining consciousness and more organized resistance may display difficul- 
ties for occupants. For long times, the resistance probability decreases more slowly. 
This means that resistance against occupation extends to long times in practice. At 
this point, the number of soldiers is also important, because the density of soldiers 
determines the speed of decaying. 

Of course, the mechanism considered in this work can be regarded as simple, but, 
it would be useful to understand the time evolution of the resistance probability of the 
community against to occupation in the one-dimensional model under some consider- 
ations. We noted that more realistic sociophysics model for better understanding the 
time evolution of resistance against to occupation appearing may be devised on the 
randomly diluted square lattice, where every site is either empty or carrier one spin 
which represent a host people or soldier randomly. Also, in that model both host people 
and soldiers move on the lattice via diffusion. Hereafter, also the sociophysics model 
mentioned above will be treated in the near future. 

Finally, it is surprisingly seen that there is a connection simple sociophysics model 
of which details are given in this paper with one-dimensional random walk processes 
in the presence of a periodic distribution of traps |18l 1191 OU] . In the both models, the 
resistance probability for any host people and the survival probability for any walker 
are given by Eq.Q at any time. Since the time evolution of probabilities in these 
models arise from the same result, hence, one can suggest that the sociophysics model 
herein and random walk model on the trapping lattice are likely to be in the same 
universality class. 
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